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Abstract 

Some power series representations of the modified Bessel functions 
(McDonald functions K^) are derived using the little known formal- 
ism of fractional derivatives. The resulting summation formulae are 
believed to be new. 

1 Fractional derivatives 

There are several non-trivial examples in mathematics when some quantity, 
originally defined as integer, can radically extend its original range and as- 
sume fractional or even continuous values. The most common example is 
the gamma-function of Euler which is a natural, unique generalization of the 
ordinary factorial: 

n! = n z = r (n + 1) = J re-'dt {n > -1) 

i=l Q 

The same thing may be performed with the order of derivatives which can 
also be made fractional. Although useful, fractional derivatives do not how- 
ever create any essentially new calculus being rather some sort of particular, 
relatively simple, integral transforms. 
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Following Oldham and Spanier (see e.g. [Q]; cf. also [Q, and [Q) we 
define the fractional derivative d^_^ = (^^j by an integral representation 
known as the Riemann-Liouville integral. Given a real number s < 0, define 

dLaf (X) = ^773^ fix- t) — ' f (t) dt (1) 

a 

where a < x is a fixed number, referred to as the boundary point. Since 
s < 0, the integral is convergent, provided that / behaves well. For s > we 
define 

where n is a positive integer chosen large enough so that s — ?t. < in order to 
assure convergence of the integral in definition It is not difficult to show 
that, as expected, such a definition does not depend on n. One can further 
prove that the familiar Leibniz rule for product differentiation has the form 

d:.aU9)=i:i']dl-if ^x^ag (2) 
j=0 

where the generalized binomial is 

fs\ ^ r(i + .) ^ (-i)^r(j-.) 

y j\T{l + s-3) j\ T{-s) 

It can also be shown that the following rule for the arbitrary order derivative 
of the power function holds 

d:.a - = ^^^l^^^ , (x - af-^ (x - a)^"" (3) 

T [p + 1 — s) [p — ny. 

for any s and p > —1. However, the analogous formulae in the case of 
the exponential and logarithm functions may, on the first sight, be quite 
surprising 

„„„ f PS / 7 {-S, Px) s=n nn 



dx' 



exp (px) = exp {i3x) / , ^ /3" exp {(3x) 



, x~ 



Inx 



dx' T(l-s) 



In X — -0 (— s) — C H — 
s. 



^ x" 



(Here n denotes a positive integer. For brevity we have set the boundary 
point a = 0.) 
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2 Some definite integrals as fractional deriva- 
tives 

Let us now consider the following definite integrals (cf. |p, formulae 3.471.4 
and 3.471.8) 



X ~ tY ^ exp 



t / a/ttx V 2x . 



(8) 

which are valid for x > 0,Re/3 > 0,Re/i > 0. On using (|I]) it is evident that 
both (P and (j^) may simply be interpreted as fractional derivatives of the 
appropriate functions. Indeed 



dx' 

QS 



V X 



dx^ 



X 2 exp 



/9 



X 



(6) 

(7) 



It is thus possible even to define the functions Ky as fractional derivatives 
with zero boundary point. 



3 Polynomials Vk 

In order to evaluate the left hand side of (|^) and let us introduce 



{(3x-) ^ x' exp ^ exp 



(8) 



It is easy to prove that, contrary to their appearance, Vk are just polynomials 

k 



j=0 
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with coefficients given by 



^ ' &^!(J-^)! r(-m) 
In particular, for a = — 1 we have 

,(-1) ^ . .^.f (-1)' r(fc + ^) ^ (fc-i)! 

^ ^ i^o^'(j-^)! r(z) (A:-j)!j!(j-l)! 

(The last equality may be rigorously proved using elementary methods pre- 
sented in e.g. 0; Mathematica 3.0 effectively simplifies such sums.) We shall 
also need fractional derivatives of the following expressions 

x''exp(-/?x") 

with a, V real. Using Leibniz rule (||) and the property (0) we have 



& [^^ exp ^ .-^'-^ exp (-,.^) E^ ./j^;:^,) ^-;"^ 

(9) 



4 Series expansion of the Modified Bessel func- 
tions 

Inserting @ into (P) or (|^) we can get a general expansion of the modified 
Bessel functions (called also McDonald functions; cf. e.g. 0) holding for 
any positive a 



KJz) 



{2z 



. , r(2.) ^(-i)^r(fc + i-. 

exp i-z) — ^ 



Vt'\2z) (10) 



Substituting polynomials (H) into ( |10|) we further get, after some manipula- 
tions 



Ks {z) = (s) Z-' exp {-z) 



oo (!_,) 



E 



'A; - 1^ 



-2z)^ 



where 

T{k + a) 



Via) 



denotes the Pochhammer symboL The same method apphed to (|^) gives 



(11) 

which is yet another power series convergent to the McDonald function. The 
famihar ratio of Pochhammer symbols in (|I0|) and (|lT]) allows us to call these 
expansions 'hypergeometric-like' since they formally resemble the well-known 
hypergeometric Kummer function iFi. 
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